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1 INTRODUCTION 



ABSTRACT 

We predict the imprint of linear bubbly perturbations on the polarization and tem- 
perature anisotropies of the cosmic microwave background (CMB). 

We analytically model a bubbly density perturbation at the beginning of the ra- 
diation dominated era and we apply the linear theory of cosmological perturbations 
to compute its time evolution. At decoupling, it uniquely signs the CMB polarization 
and temperature anisotropy sky. As predicted by recent general work regarding spa- 
tially limited cosmological seeds, during evolution the perturbation propagates beyond 
the size of the bubble and reaches the CMB sound horizon at the time considered. 
Therefore, its signal appears as a series of concentric rings, each characterized by its 
own amplitude and sign, on the scale of the sound horizon at decoupling (< 1° on 
the sky). Polarization and temperature rings are strictly correlated; photons coming 
from the center of the bubble are not polarized, because of the spherical symmetry 
of the present problem. As expected for linear perturbations with size L and density 
contrast 6 at decoupling, 8T/T is roughly 6(L/H ) 2 ; the polarization is about 10% 
of the temperature anisotropy. 

We predict the impact of a distribution of bubbles on the CMB polarization and 
temperature power spectra. Considering models containing both CDM Gaussian and 
bubbly non-Gaussian fluctuations, we simulate and analyze 10° x 10° sky patches 
with angular resolution of about 3.5 . The CMB power associated with the bubbles is 
entirely on sub-degree angular scales (200 < I < 1000), that will be explored by the 
forthcoming high resolution CMB experiments with the percent precision. Depending 
on the parameters of the bubbly distribution we find extra-power with respect to the 
ordinary CDM Gaussian fluctuations; we infer simple analytical scalings of the power 
induced by bubbly perturbations and we constrain our parameters with the existing 
data. 



et al. 1998 and references therein); they shall culminate with 



World wide theoretical and experimental efforts to under- 
stand what we could learn from the Cosmic Microwave Back- 
ground I (CMB) ate presently in progress. Anisotropies in 
the CMB shuuld euiilaiii pieuiuus details abuuL high en- 
ergy physics processes occurred in the early universe, still 
hidden to our knowledge. Reasonably, these processes left 
traces that have been stretched out to large and observ- 
able scales by an early era of accelerate expansion, called 
inflation; at decoupling between matter and radiation, these 
traces imprinted anisotropies in the CMB, thus being possi- 
bly detected and recognized. 



the Planck mission of the next decade, that will provide the 
whole sky anisotropy map down to a minimum ST/T ~ 10 -6 



and an angular resolution of about 10 ( jee the Planck Sur- 



veyor and Microwave Anisotropy Probe home pages) 



According to the simplest inflationary phenomenology, 
at very high energy in the early universe (T ~ 10~ 5 Planck 
masses), a scalar field slowly rolls toward the minimum of 
its potential, giving the non-zero vacuum energy responsible 
for the expansion itself. The quantum fluctuations gener- 
ated during this phase are stretched out by the inflationary 
expansion itself to large scales and become the perturba- 



Recently, the treatment of the linear CMB inhomo- 


tions we observe today in the matter power spectrum ( 


set 


geneitk 


s has been put in a complete and organic form (|Hu| 


Mukhanov, Feldman & Branderberger 1992 for an extens 


ivc 


et al. 1 


997|). At the .same time, many experiments arc at 


treatment!). However, this is not the only known inflationary 



work to explore the C MB anisotropies toward smaller and 
smaller angular scales ( see Lasenby et al. 1998, De Bernardis 



mechanism. Many fundamental fields may act on stage and 
the effective potential may have several minima separated 
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by potential barriers. If this is the case quantum tunneling 
occurs, and the nucleated bubbles are stretched out to large 
scales just like the ordinary quantum fluctuations. At reheat- 
ing the energy stored in the bubble shell is converted into 
matter and radiation and a bubbly trace is left in the den- 



Starting from the signal from isolated bubbles, we pre- 
dict the effect of a distribution of bubbles on the CMB 
temperature and polarization power spectra; as expected 
in scenarios admitting a bubble nucleating era during the 



context of first order inflation ( see Kolb 1991 and references 
therein^; the possibility that bubbles have left traces in the 
large sc ale matter distribution was considered in the las t 
decade ( ^ee Amendola et al. 1996 and references therein ), 
in connection with the increasing obser yational evidence of 
large v oids in the galaxy distribution (El- Ad, Piran & da 
Costa, |L997| ). 

The aim of this work is to accurately predict the CMB 
polarization and temperature anisotropies from bubbly den- 
sity perturbations existing at decoupling, relics of the infla- 
tionary tunneling phenomena. In an early work ( Baccigalupi 
et al. IS 197) we considered anisotropies from non- linear bub- 
bles mainly for what c oncerns metric in duced perturbations. 
Then, more recently (Baccigalupi 1998) we used an approx- 
imated method to compute the temperature anisotropies 
ST/T from isolated inflationary bubbles in linear regime; 
we found that, just like a pebble in a pond, the bubbly 
perturbation generates anisotropy waves, traveling outward 
with the CMB sound velocity and reaching the scale of the 
sound horizon at the time we are examining it. Therefore, 
the image of an inflationary bubble at decoupling is a se- 
ries of concentric isothermal rings of different color (sign of 
ST/T) on the scale of the sound h orizon at decoupling (< 1° 
in the sky). In a separate work (Amendola, Baccigalupi & 
Occhionero 1998) we explored the impact of a distribution of 



inflationary s low roll (see Amendola et al. 1996 and refer- 



sity dis tribution. Tunneling processes are predicted in the ences therein ), we consider models containing a population 



bubbles on the CMB temperature angular power spectrum, 
finding for a wide set of statistical parameters extra power 
on sub-degree angular scales. 

Here we compute the CMB polarization anisotropies 
induced by inflationary bubbles, and we generalize previous 
results on the temperature in the cont ext of the full linear 
theo ry of cosmological perturbations (Ma & Bertschinger 



1995|^JVe model the remnant of an inflationary bubble as 
a linear bubbly perturbation in the density distribution at 
the beginning of the radiation dominated era. Each Fourier 
mode is analytically calculated and numerically evolved fol- 
lowing the linearized Einstein and Boltzmann equations. We 
use these results to study the signal from a single bubble. 
Regarding the size and density contrast of the bubbles, they 
are ultimately determined by the epoch of nucleation dur- 
ing inflation and by the properties of the tunneling barrier, 
so that these para meters depend on the details of the infla- 
tionary potential ( see Amendola et al. 1996 and references 
therein 
Baccig; 



as in previous works ([Baccigalupi 1998 



Amendola. 



properties correspon ding to the large voids observed in the 
galaxy distribution (El- Ad, Piran & da Costa, 1997), that 
have deep cavities Sp/p ~ 1 and comoving radii of about 
20^ 40ft" 1 Mpc; if they are of primordial origin, the evolu- 
tion from decoupling to the present implies an overcomoving 
growth by a factor of 2, and the density cont rast increases 
approximatively linearly with the scale factor ( Dcchionero et 
al. 199^) : thus, the voids corresponding to the present ones 
should have at decoupling a comoving radius R = 10-f-20/i -1 
Mpc and density contrast S = 10 -3 -r- 10~ 2 , values that we 
adopt in this work. 



of bubbles together with ordinary CDM Gaussian perturba- 
tions; in this context we simulate CMB sky maps and we 
compute the anisotropy angular power spectrum for various 
bubbly distribution parameters. As a general useful com- 
ment here, we wish to emphasize that both in the cases 
of a single bubble and a distribution of them, the signals 
treated here are definitely different from Gaussian pertur- 
bations. Gaussianity means that the fluctuations exhibit a 
completely disordered spatial distribution; only in this case 
the CMB angular power spectra tell us everything about the 
signal. On the contrary a bubbly perturbation is an ordered 
structure: it really possesses a spherical symmetry. Techni- 
cally speaking, its signal is non-Gaussian; therefore, even if 
the CMB power spectrum is a powerful tool to investigate 
also non-Gaussian models, it does not fix completely the un- 
derlying fluctuations, and other methods should be adopted 
to investigate which type of non-Gaussian ordered pattern is 
present on the observed map. The general treatment of CMB 
anisotropies from cosmological structures having spatial fi- 
nite ex tension and symm etries has been recently treated in 
detail (Baccigalupi 1999); we apply the general treatment 
developed in that work to the present case. 

The paper is organized as follows: in section II we recall 
the relevant equations for our purposes and we expose the 
computational details; in section III we evolve the bubbly 
perturbation to study the corresponding CMB perturbation 
at different times and to predict the effect on the CMB polar- 
ization and temperature anisotropies from isolated bubbles; 
in section IV we consider several distributions of bubbles and 
we predict their imprint on the CMB angular power spec- 
tra, also giving simple scalings in terms of the distribution 
parameters; finally, section V contains the conclusions. 



2 THE EQUATIONS SYSTEM 

The evolution equations for fluid and CMB quantities may 
be obt ained by the Boltzmann a nd linearized E instein equa- 



tions (Ma & Bertschinger 1995; Hu et al. 1997). A standard 
CDM flat background is assumed, including cold dark mat- 
ter ( c ), baryons (&), photons ( 7 ) and three families of mass- 
less neutrinos („). The background flat Friedmann Robert- 
son Walker (FRW) metric is 



lupi & Occhionero 1996) here we choose the bubble ds 2 = a{rf} 2 ( 



~dr] 2 + dr 2 + r 2 dfi 2 ) 



(1) 



where rj(t) — ^dr/a(r) is the conformal time. The back- 
ground evolution is driven by the Einstein equation 



87rG 2 \ ~* 



(2) 



where the index x runs over all the fluid species; the back- 
ground parameters are fio = 1, h = .5, fit = 0.05, Qcdm = 
1 — fib. The perturbed metric tensor is 



(3) 
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where a(r;) 2 7 MI , represents the background. Since /i M „ <C 
7 M „, a gauge freedom reduces the number of physically sig- 
nificant quantities in the perturbation metric tensor; in this 
work we adopt the generalized Newtonian gauge in which 
the two scalar perturbed m etric component are = Hqq/2 
and 2$ = hn = /122 = ^33 (Ma & Bertschinger 1995); con- 
cerning the fluid, 5 X indicates the density contrast 5p/p for 
the species x and v x its peculiar velocity in the background 
metric. Unless otherwise specified, all the equations for the 
perturbed quantities in the remaining part of this section 
are written in Fourier space, and the wavenumber argument 
is omitted. 

The equations for the matter species are: 



-kv c - 3fc $ 



-kvb—3k $ , i>t 



k^ 



a , . 4p 7 . . 

-— U6+KW + - — an e aT(v-y—v b ) 
a ipb 



(4) 



,(5) 



where f = ax e n e aT is the differential optical depth; gt is the 
Thomson scattering cross section, n e is the electron number 



density and x e the ionization f raction (see Hu & Sugiyama 
1995 fo 



useful fitting formulas) 



Now we give the expressions for the CMB temperature 
and polarization perturbations in the context of a spheri- 
cal pertu rbation field; thes e formulas have been recently de- 
veloped (Baccigalupi 1999). As customary, we indicate the 
temperature perturbation 5T/T with O and the linear po- 
larization amplitude with the Stokes parameters Q, U. 

Temperature perturbation carried by CMB photons 
scattered on an arbitrary direction h at a spacetime posi- 
tion (rj, r) around a spherical structure is given by 



e 



5>c 



k 2 dk 
2tt 2 



(6) 



where Pi and ji are respectively Legendre polynomials and 
fractional order Bessel functions. All the specifications re- 
garding the perturbation source are encoded in the Qi coef- 
ficients, which obey evolution equations that we shall write 
in a moment. Concerning the polarization, we have to spec- 
ify the frame axes describing the polarization tensor on the 
plane perpendicular to n. We choose them parallel and per- 
pendicular to the plane formed by ri and f; with this choice, 
the pola rization tensor is described only by Q ( Baccigalupi 
1999) iJi the following way: 



E 



(l- 2)! 
(1 + 2)1 



Pi\n- 



k 2 dk 
2tt 2 



E l ( V ,k)j l (kr) 



(7) 



Pi are second order Legendre polynomials; this is a sig- 
nificant difference from (^): correctly, they make photons 
moving radially not polarized (P 2 (±l) = 0). Again all the 
characteristics of the perturbation source are contained in 
the Ei quantities. Now we give the evolution equations for 
the Ei and <di coefficients. They are the amplitudes of po- 
larization and temperature CMB pertu rbations expand ed in 
scalar and tensor spherical harmonics (Hu et al. 1997), and 
obey the following equations: 



e = --e 1 



6i = fc6 - -fc6 2 + r(v b - 6i) + h% 



(8) 



(9) 



E 2 = -^kE 3 -f(±e 2 + lE 2 ] 



and for I > 3 
I 



@i = k 



Ei = k 



21-1 



©1-1 - i^4©i+i 



VP 



21-1 



21 + 3 

V(l + l) 2 -4 
21 + 3 



E1+1 



+Ei 



(10) 
(11) 

(12) 
(13) 



In Newtonian gauge the lowest multipoles are linked to 
the photon fluid quantities by S 1 = 40o, f 7 = 0i and 
7r 7 = 126>2/5. Massless neutrinos can be treated as photons 
without the polarization and Thomson scattering terms. 
There remain the equations for the gravitational potentials: 



4nGa (p c 8 c + PbSt + P-ySj + p v 5 v ) + 
2 3 a 

i n r.. hi n > 1 * 



+ 4nGa 



3d/ 44 
k a \ pcVc + PbVb + oP~t v ~tvP vVv 



, 2 / -r t \ 8irG . . 
A; (* + $) = — — (p 7 vr 7 + p„ir u ) 



(14) 



(15) 



As it is known (Ma & Bertschinger 1995), at early times 
the above system can be solved by using the tight coupling 
approximation between photons and baryons. The multipole 
equations are expanded in powers of klf <C 1. The only 
zero order terms are Oo and &i from (plpl|9|), and obey the 
following equations: 

(16) 



9 
drj 



3pt 
4p 7 



9i 



fc9 + k 1 + 



3pb 
4p 7 



(17) 



where Oi is assumed to coincide with Vb to the lowest order. 
Increasing the order in k/f the higher order multipoles are 
given by 



e 2 



k 8 
?9 



6i 



/ 



-e. 



E-2 



Ei 



H 2 



k VP 



-Ei. 



(18) 
(19) 



Ql i2l~l" l ~ L ' i 21-1 

We integrate in time the system (^Jl^ @Jl|) until k/f = 
.1 occurs, thereafter integrating the complete equations; of 
course, care is taken that the results do not depend at all 
on this choice. 

We take adiabatic initial conditions: at early times 
5 C — 5b — 3<5 7 /4 = 3<5„/4 (all the velocity are initially zero) 
and the second member in equation (n4) at r\ — is pro- 
portional for each Fourier mode to the initial perturbation 
spectrum, that we define now. As we mentioned in the in- 
troduction, a well understood property of nucleated bubbles 
in field theory is their spatial growth; once nucleated in a 
de Sitter background like the inflationary one, a bubble is 
stretched out to super-horizon scales like the ordinary slow- 
rolling Gaussian fluctuations. Thus, although the coupling 
between the inflationary fields and the known sector of parti- 
cle physics is unknown, a reasonable and simple assumption 
is that at the end of inflation bubbles simply convert their 
energy density distribution in mat ter and radiation, leaving 
a bubbly trace in the cosmic fluid ( see Amendola et al. 1996 
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and references therein) . We parameterize the latter with the 
f ollowing model, th at we already adopted in an earlier work 
( |Baccigahrpi 1998[ ): 

p \ ~ 2, 

Sp 
P 



— = A + B sin 



2tt . 

— [x + a 
a 



1) 



2<,<1 



(20) 



S -P=C + Dsin 
P 



2tt , 

(X + G 

a 



1 



1 < x < 1 + - 



x = r/i? is an adimensional radial coordinate with J? being 
the comoving radius of the bubble, and a is the shell thick- 



The Fourier transform (Sp/p)k = J (dp/ p)(r)t 



ik-r i3 



d 6 x = 



(i )(sin h /ki )4;i i di — uf Ik 



ui L m baLiun 



deep oscillation that appears at the indicated redshift; also, 
a positive crest arises at a distance corresponding roughly 
to the radius of the bubble. At the time indicated in panel 
(c) an opposite oscillation occurs in the center, forming a 
negative crest, while the positive one is traveling outward. 
Finally, in panel (d) we can see the situation just before de- 
coupling, thus resembling the real CMB anisotropics that 
we shall show in the following. A central negative spot is 
left, together with a well visible wave, with negative and 
positive crests, that is traveling outward with the sound ve- 
locity and that will be photographed by the CMB decou- 
pling photons. Up to an accuracy of percent, the relevant 
contribution to the perturbation is given by the lowest mul- 
tipoles Oo and Qi - that we have us e d for the computations 
in previous works (tBaccigalupi 1998j ; Amendola, Baccigalupi 



easily c alculated analytically. The apcctrum ia of courac 

strongly scale dependent: all the power is essentially concen- 
trated around the wavenumber corresponding to the bubble 
radius k ~ 2n/R, while oscillations at higher k describe the 
bubble shell. The constants are determined by the requests 
of continuity and compensation: 



D 



C = B 



A + 



(2-a) 2 + 8- 



2 

2(4- 



a)/n 2 



(21) 



16 + 4cr 2 - 24a 3 / 71 " 2 
The physically most relevant parameters are 8 and R: the 
former sets the amplitude of the perturbation and the latter 
its comoving size. The shell thickness does not change the 
substance of the result, and from now on we fix a — .3 . 
In order to make the following results more clear, S is nor- 
malized with the density contrast taken in the center of the 
bubble at decoupling between matter and radiation. 

The very final thing to say is about the initial condi- 
tions in the CMB equations (|||| [u|l^,|l|Jl3|) ; everything is 
initially zero except for the lowest multipole of the temper- 
ature perturbation (Ma & Bertschinger 1995) 



8o(0) = -2* (0) 



(22) 



thus being proportional, for each Fourier mode, to the initial 
perturbation spectrum. 

In the next section we integrate the equations defined 
here to see how the bubbly perturbation evolves in time and 
how it would appear in an high resolution CMB anisotropy 
observation. 



3 TEMPERATURE AND POLARIZATION 
EFFECTS FROM A SINGLE BUBBLE 

Let's put the bubbly perturbation (^o|) as initial condition 
of the equations and let it evolve. We are interested in the 
behavior of the CMB perturbation, so we take some photos 
at different times. First, look at figure (1). It represents the 
pure temperature perturbation induced by a bubble with 
radius 7? = 20/i _1 Mpc at the indicated times expressed in 
redshift; on the z-axis we can read the comoving distance 
from the center, and the photon propagation direction is or- 
thogonal to the radial direction, n-f = 0, as indicated. Panel 
(a) shows the initial condition, that remains unchanged until 
the horizon crossing, that for the scale that we have consid- 
ered occurs nearly at equivalence. Panel (b) shows a central 



fc Occhionero 1998D; in each panel we have plotted the con- 



tributions from these two terms only (thin dashed lines); 
the difference is not visible since they agree at the level of 
accuracy specified. 

Remarkably, the same phenomenology occurs for polar- 
ization. In figure (2) we plot the polarization amplitude Q 
defined in the previous section as a function of the comoving 
radial distance r and for propagation perpendicular to the 
radius. The initial condition (panel (a)) is vanishing since 
the first multipole of the polarization components is I — 2 
and it is coupled with the same momentum of the pure tem- 
perature perturbations, that is zero at the beginning; phys- 
ically, the reason is that polarization is coupled to acoustic 
oscillations that do not occur when all the modes are out- 
side the effective horizon. At the horizon crossing, panel (6), 
a negative oscillation occurs, followed by an opposite one 
(panel (c)); note that the light coming from the center is 
not polarized, simply because light propagating radially in 
a spherically symmetric problem has to be unpolarized for 
a pure geometric reason. Just before decoupling, a polariza- 
tion wave is well visible (panel (d)); it is placed at the same 
distance from the center as the Q wave in figure (1); this is 
an expected feature, due to the coupling between the polar- 
ization and pure temperature perturbations, manifest in the 
system §|Jlol0, 

The general phenomenology is therefore similar to the 
waves in a pond when a pebble is thrown into it. This is a 
general behavior valid for any c osmological struc ture hav- 
ing a spatially limited extension (Baccigalupi 1999). We can 
physically understand this point by considering the equation 
for the / = multipole of the temperature perturbation (t o 
the lowest order in k/f). As it is known (Hu et al. 1997), 
from (|16[|l7|) we get 



drj \4 Pl 



3r : 



d f 3p b 
drj 



4> 



(23) 



where c s = 1/ + 3pt/4p 7 ) is the sound velocity of the 
photon-baryon fluid. The right hand side contains the grav- 
itational forcing terms, while the left hand side is just a 
damped wave equation. For a localized initial inhomogene- 
ity, the corresponding CMB perturbation propagates beyond 
the initial size, generating waves traveling outward with the 
sound velocity, and reaching the size of the sound horizon 
at the time we are examining it. 
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Let us consider now the numbers. The central den- 
sity contrast at decoupling 8 has been factored in the y- 
axis labels of the figures. Roughly, the amplitude of 
follows the known esti mates existing for linear perturba - 
tions, Q ~ S(R/ H^ 1 ) 2 (see Padmanabhan 1993 for reviews), 
where the Hubble radius is calculated at decoupling. The 
polarization amplitude is roughly 10% of the temperature 
perturbation; this is also an expected result, since the lowest 
multipole of the polarization is coupled with the quadrupole 
momentum of the temperature perturbation (O2) that in 
turn is a first order effect in k/f with respect to the quan- 
tities in ( p^|Jl7| ), as is evident from (|l8|); the same consider- 



ation holds for the ordinary Gaussian perturbations (Hu et 
al. 199^). 

The graphs we have shown until now were merely photos 
of the bubbly CMB perturbation at different times. Now we 
want to predict in detail the CMB anisotropies. Unlike Gaus- 
sian scale-invariant perturbations, our bubbly CMB signal 
is localized and vanishes beyond a sound horizon from its 
center. Therefore, its position along the photons path is an 
adjunctive degree of freedom with respect to the ordinary 
perturbation models and it has to be specified. Because of 
the spherical symmetry of the bubble, it can be simply the 
distance d between its center and the last scattering surface 
(LSS) peak; the latter is defined as follows. CMB photons 
have probability P(rj) — fe~ T (r is the differential optical 
depth defined in the previous section) to be last scattered 
between z and z + dz; the detailed form of P{ri) depends on 
the particular cosmological scenario of course, nevertheless 
it is in most case s a Gaussian centered around z ~ 1100 with 
width Az ~ 50 (see Hu & Sugiyama 1995 for useful fitting 



formula]!). We define the LSS peak as the point from which 
CMB photons were last scattered with highest probability 
(maximum P(ri)) along the direction corresponding to the 
bubble center in the sky. Once d is specified, we can get the 
CMB anisotropies via simple line of sight integrations: 



Q(now, here, h) = 



Q(now, here, h) 



(9 + ^)[r ) ,r{ n ,h),h]P{r l )dr ] , (24) 



'\q, r(r), n), n]P(rf)drj 



(25) 



where the Sachs-Wolfe effect on photons climbing out of the 
potential hills g enerated by the bubble, simply given by *P 
(Hu et al. 1997), has been added; d, n and the propagation 



direction n fix the radius r in which we have to compute the 
CMB quantities in the integrals: as it is easy to see 



r{ri,h) = [(d + r]a — tjlss) 2 + 



+ (Vo - v) 



2(d + r/o - riLss)(-no - v) n ■ n, 



,1/2 



(26) 



where n c is the direction of the radiation coming from the 
center of the bubble. In fact, since rjo — r\ is just the causal 
distance covered by a photon last scattered at rj and reaching 
us today, equation (Eq) is obtained with simple trigonometry 
(Baccigalupi 1999). Also we can define naturally the useful 
angle 9 by 



n ■ Tic = cos ( 



(27) 



it is simply the angle between the photon propagation direc- 
tion h and the direction corresponding to photons coming 
from the center of the bubble in the sky. 

Now we are ready to see the results. Figure (3) shows 



the CMB polarization (top) and temperature (bottom) 
anisotropies for a bubble at decoupling and with the size 
indicated; signals are plotted as a function of 6. 

Look at the top panels first. They show the polariza- 
tion anisotropy Q of the light coming from the bubble; we 
recall that it represents the difference between the temper- 
ature fluctuations polarized along the axes parallel and per- 
pendicular to the plane formed by the radial and the pho- 
ton propagation direction. At small 9, Q is vanishing as it 
must be, since the radiation coming from the bubble cen- 
ter is propagating radially and thus must be not polarized. 
The polarization wave that we have shown in figure (2) has 
been photographed by the decoupling photons here: it is well 
visible for both the bubble sizes chosen at an angular posi- 
tion corresponding to the sound horizon at decoupling. Note 
that the bubble is confined to small angles, say at about 
9 < R/2Hq 1 ~ 10 for the radii treated here. Therefore, 
the polarization wave beyond the bubble size is a unique 
sign of its presence. Also note that in each panel we have 
plotted three curves. The solid one correspond to a bubble 
exactly lying on the LSS peak (d = 0), while the short and 
long dashed ones correspond to a bubble lying backward 
(d — —R) and forward (d — +R) with respect to the LSS 
peak: these curves allow us to see the dependence of the 
anisotropy on the relative disposition of the perturbation 
source with respect to the LSS; of course, if d was greater 
than the sound horizon, no signal at all could be visible. 

Look now at the bottom panels. It is as before, except 
for the fact that here is plotted the temperature anisotropy. 
Note the central negative spot, at the same location of 
the bubble; again, at a sound horizon from the center, the 
ST /T wave propagating outward with the photon-baryon 
fluid sound velocity has been photographed by the last scat- 
tered p hotons. We have previously computed the curves for 
ST/T ( Baccigalupi 1998 ) by considering only the monopole 
and the dipole (Oo and 0i) contributions to the signal (thin 
continue line); here we obtain again the previous results by 
considering all the multipoles, showing at the same time 
their good accuracy. 

For what concerns the numbers, a comparison with fig- 
ures (2) and (3) shows how the CMB anisotropies are weaker 
than in the images of the perturbation at different times. To- 
gether with the partial cancelation of the opposite signs of 
O and 9, this is also because of the line of sight integration 
across the LSS; the size of the bubbles is comparable with 
the LSS width, and due to the oscillating behavior of the 
signal, the effect of the integral is a certain weakening of the 
resulting signal. 

Now, anticipating issues of the next section, we realis- 
tically predict what we should observe in the CMB sky if 
a bubble was really present at decoupling on the observed 
patch. In figure (4) we have simulated 2° x 2° portions of mi- 
crowave sky. First look at the bottom panel. In each point we 
have plotted the Q Stokes parameter of the CMB anisotropy. 
The signal consists of two compon ents: a Gaussian scale in - 
variant adiabatic CDM spectrum ( |Ma fc Bertschinger 199 



and a bubble with comoving size R — 20h~ L Mpc lying ex- 
actly on the LSS (d — 0), just like the solid line of figure 
(3). This case could be realistic in the sense that if during 
inflation more than one field were on the scene, and at least 
one undergoes quantum tunneling and bubble nucleation, at 
reheating we should expect both the kinds of perturbation 



© 0000 RAS, MNRAS 000, 000-000 



6 C.Baccigalupi & F.Perrotta 



of the figure, bubbly and Gaussian; also, they simply add the sky via the Heal pix pixelization (see the home page oi 



together if the perturbations are linear. The fascinating cir- 



the Healpix package), setting the signal outside the patch 



cular imprint of the bubble is clearly visible well over the 
CDM contribution; the central density contrast in this case 
is 5 = 10~ 2 . Note that the bubble lies on the center of the 
figure, within 10 from the center. The circular polarization 
wave is well visible and it is outside the bubble, as we have 
already pointed out. If 8 is reduced, we go to the case of the 
upper panel, in which the bubble is also present, but with 
a smaller amplitude, 8 — 5 • 1CP 3 . In this unlucky case the 
detection of the bubbly signal would require the use of more 
sophisticated image analysis tools than the human eye. 

Finally, figure (5) shows the temperature signal instead 
of the polarization. Note as the circular rings are well ev- 
ident in the bottom panel, where 8 is higher. Again, the 
bubble lies in the very central region, corresponding to the 
small central dark spot; the rings around are a feature of the 
CMB own physics and are coupled, for extension and oscil- 
lating behavior, with the polarization ones. The central dark 
spot is absent in the polarization case due to the symmetry 
constraints previously mentioned. 

The cross detection, both in temperature and polar- 
ization anisotropics, of the CMB rings caused by a local- 
ized perturbation like a bubble would be a powerful check 
if signals like the present ones should be really detecte d in 
the fut ure high resolution maps from MAP and Planck ( sec 



the Pla lck Surveyor and Microwave Anisotropy Probe home 



pages^T 

In the next section we consider the imprint of a popu- 
lation of bubbles on the CMB anisotropy power spectrum. 



4 TEMPERATURE AND POLARIZATION 
EFFECTS FROM A DISTRIBUTION OF 
BUBBLES 

Let us consider now a distribution of bubbles. Our purpose is 
to check whether, and under which hypothesis, a distribution 
of bubbles is able to make a distinctive imprint on the CMB 
angular power spectrum: let us spend a few words to recall 
its definition. Given a sky signal s(9,<f>), the angular power 
spectrum is a set of coefficients Cis simply defined in terms of 
the coefficients a; m of its expansion into spherical harmonics: 



Ci = 



m — l 

1 V- | ,2 

dim 



21 + 1 ^ 



(28) 



-l 



Since it depends on the squared modules of the ai m s coeffi 
cients, it is sensitive on the anisotropy power on th e angular 



scale ol roughly correspond ing to 180/Z degrees (see Pad- 



manabhan 1993 for reviews 



This work contains some improvements with respect to 
previous results treating the CMB power spectrum from 



a distribution of bubbles (Amendola, Baccigalupi & Oc- 



chionero 1998): first both polarization and temperature sig- 



nals are computed here, and second the signals from far 
bubbles, namely up to a sound horizon from the LSS, are 
considered instead of only the closest ones at a distance of 
about a radius from the LSS. 

We analyze 10° x 10° zero mean CMB anisotropy maps 
with 3.5 x 3.5 squared pixels. Bubbly signals are ran- 
domly distributed within the patch, which is projected on 



to zero; successively, the signal is developed into spherical 
harmonics using Fast Fourier Transform methods in order 
to get the Ci coefficients. We analyze both maps containing 
only the signals from bubbles and the latter together with 
the CDM Gaussian perturbations as in figures 

First, in order to test our analysis method, we consider 
only CDM fluctuations. Starting from the theoretical CMB 
power spectrum, we simulate a patch as specified above and 
we extract from the latter the power spectrum itself, thus 
reconstructing the theoretical one. In figure ^ we plot the re- 
constructed and theoretical CDM power spectra for temper- 
ature (up) and polarization (down) anisotropies. The spectra 
are in quite good agreement (roughly 10% accuracy) for a 
wide range of angular scales, corresponding to I > 50, while 
of course they disagree as approaching the dimension of the 
patch (roughly corresponding to I ~ 20); the oscillations in 
the simulated spactra are due to the sample variance on the 
sky patch, but for our purpose here this method works well. 

Let us define now the distribution of the bubbles. Be- 
cause of the anisotropy propagation treated in the previous 
section, the signal on the LSS contains contributions from 
bubbles centered at a distance that extends up to a sound 
horizon r a at decoupling. Indicating with R the radius of 
the bubbles, and with x the fraction of volume occupied, 
the number of bubbles that we should consider in our patch 
is roughly 



n = x ■ 2r s 



(2Ho 1 ) 2 dQ. 



(29) 



where dil is the solid angle subtended by the patch in the 
sky and 2H^ 1 is the comoving distance of the LSS from us. 
In our simulations we adopted R = 20/i _1 Mpc, so that for 
r s ~ lOO/i" 1 Mpc a x — 10% volume fraction is covered by 
roughly 600 bubbles in the specified patch. 

First, let us analyze the effect of the bubbles only. Fig- 
ure shows the temperature (up) and polarization (down) 
spectra from bubbles filling a volume from 10% (thin line) 
to 50% (heavy line). The bubbles have central density con- 
trast at decoupling 8 ~ 10~ 3 . As expected for sub-horizon 
structures at decoupling, the whole signal lies on sub-degree 
angular scales, I > 200. In the temperature case, two re- 
gion can be identified: the highest one describes the central 
hot spot of the bubbly signals and involves the correspond- 
ing multipoles at I ~ 800; while the signals vanishes rapidly 
higher is, it shows a sort of tail that reaches I ~ 200, and 
that is generated by the anisotropy sound waves treated in 
the previous section extending the signal nearly up to one 
degree, see figures The spectrum form is completely 

different in the case of the polarization. This is an expected 
feature, since the polarization signal described in the previ- 
ous section does not possess the central hot spot due to the 
spherical symmetry of the bubble. Instead, the polarization 
spectra contains peaks roughly at four angular scales that 
are related to the polarization waves departing from the cen- 
ter of the bubbles and extending up to about one degree and 
that are evident in figures 

Let us come now to the evaluation of the scalings of the 
CMB power induced by bubbles. In this case all the bub- 
bles have the same central density contrast 8, so that their 
power spectra scale exactly as 8 2 . Moreover, as expected for 
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spatially uncorrelated bubbles, the amplitude is quite linear 
into the occupied volume fraction x. In this work we do not 
show graphs relative to other radii, but si nce the bubbles 



are su b-horizon struc tures at decoupling (see Padmanab- 



han 191)3 for reviews), we checked that their signal scales 
roughly as R ; in addition, as we mentioned in the previous 
section, the computation of the anisotropies causes an ad- 
ditional damping due to the line of sight integration across 
the LSS: the signal therefore drops by a factor 10 passing 
from 20ft -1 Mpc to 10/i -1 Mpc, implying a scaling of about 
R 3 5 . Putting together all these ingredients, and defining the 



bubbly spectral power as C 



Z(Z + l)C/'V27r calculated 



roughly at the maximum, we can infer the follow empirical 
scalings: 



C 1 



4 ■ 10"V 



R 



(50%) \2()// '.Mix- 



er = 3 • 10" 6 <5' 



R 



(50%) \2()// '.Mpc 



(30) 
(31) 



Note that these numbers, for what concerns temperature, 
lead to results that are bigge r by nearly a factor of 2 than the 



scaling inferre d previously (Amendola, Baccigalupi & Oc- 



chionerj 1998). This is probably because in our simulations 
more signal imprints the patch, since we considered the con- 
tributions from bubbles placed up to a CMB sound horizon 
away from the LSS. Since the recent measurem ents of CMB 
anisotr opy power on sub-degree angular scales ( jee Lasenby 



et al. 1998, De Bernardis et al. 1998 and references therein) 
indicate values roughly at the level 1(1 + l)C;/27r ~ 4 • 10~ i0 
in the region interesting for us, bubbles populations with 
C T larger then that values are likely to be excluded. This 
implies the following constraint on our parameters: 



c2 ( x \ I R 

V50%; I 20ft- 1 Mpc 



< 10" 



(32) 



This condition is obtained considering only the power from 
bubbles, anyway it holds also when considering the power 
from Gaussian CDM perturbations, simply the two pertur- 
bations mechanisms are uncorrelated, so that their powers 
approximatively add. Indeed, figures g and ^ show the CMB 
angular power spectra in the cases in which bubbles and 
CDM Gaussian perturbations are present in the simula ted 
patch. Also plotted are the existing data on I > 200 ( see 



Lasenb y et al. 1998, De Bernardis et al. 1998 and references 



therein] . In figure g the bubbly distribution are taken from 
the x = 50% case shown in figure M. As we anticipated, and 
as it is evident looking at figures ^jjj the power from each one 
of the two perturbations mechanism approximatively adds. 
In figure ^ as example of the effect of varying 5 is shown: 
we plotted a case in which the volume fraction is x = 10% 
but 5 = 2- 10 -2 . As a final comment, we wish to point out 
how there are cases in which the bubbles are present, filling a 
significant fraction of the volume and having properties able 
to influence the structure formation process between decou- 
pling and the present, but that produce CMB power that is 
subdominant with respect to the Gaussian CDM one: this 
happens roughly for C T < 10~ 10 , equivalent to C p < 10~ 12 . 
This is not surprising, since in these kind of non-Gaussian 
models, the power spectrum alone does not fix completely 
the underlying perturbations properties, since, as it is evi- 



dent from its definition ( |28| ) , it is sensitive on the anisotropy 
power at a given angular scale. 



5 CONCLUSION 

The most known inflationary model leaves traces in the 
form of Gaussian perturbations. With these hypothesis, it 
uniquely marks the angular power spectrum of the cosmic 
microwave background (CMB) anisotropies. However high 
energy physics may be more complicated and may leave 
other (and richer) traces, in the form of non-Gaussian per- 
turbations. During inflation, quantum tunneling phenomena 
may have occurred, and the nucleated bubbles may have 
been stretched to large and observable scales by the infla- 
tionary expansion itself. Here we predict the imprint of bub- 
bles on the CMB, both on the polarization and temperature 
anisotropies. We consider an analytical bubbly density per- 
turbation at the beginning of the radiation era, parameter- 
ized by its comoving radius R and density contrast 5. We 
perform its evolution via the linear cosmological perturba- 
tion theory. We find that it uniquely marks both the po- 
larization and temperature anisotropies. Then we consider 
a parametric distribution of bubbles and we predict its im- 
print on the CMB angular power spectra, as a function of 
the bubbles abundance and properties. 

During the evolution toward decoupling, we take some 
photos of the bubbly CMB perturbation at different times. 
As expected, we find that the initial condition remains un- 
changed until the horizon reenter occurs, and thereafter 
the perturbation starts to oscillate. This behavior generates 
sound waves propagating outward with the sound velocity of 
the photon-baryon fluid and reaching the sound horizon at 
the time we are examining the system. At decoupling, we get 
the CMB anisotropies by considering the bubbly perturba- 
tion intersecting the last scattering surface. The dependence 
of the signal on the position of the bubble along the line of 
sight has been considered. 

The CMB polarization image of an inflationary bubble 
at decoupling is a series of rings concentric with the center 
of the bubble, extending on an angular size corresponding to 
the sound horizon at decoupling (~ 1° in the sky); each ring 
is characterized by its own value of the Stokes parameters, 
giving the difference of the temperature fluctuations pro- 
jected on the axes parallel and perpendicular to the plane 
formed by the radial and photons propagation directions. As 
an important difference with respect to the case of temper- 
ature anisotropy, that reflects also in the CMB power spec- 
trum from a distribution of bubbles, the light coming from 
the center of the bubble is not polarized, since the radial 
propagation in spherical symmetry is an axial symmetric 
problem, so that no preferred axis exists on the polarization 
plane. Instead, temperature anisotropies are generally char- 
acterized by a strong central spot, together to the external 
rings. This is an important difference, that we prepare to in- 
vestigate for what concerns the impact on the CMB angular 
power spectrum. Finally, the amplitude of the CMB fluctu- 
ation follows the known estimates for linear perturbations 
with size R and density contrast 6: 5T/T ~ S ■ (R/H' 1 ) 2 
where 5 and the Hubble radius are computed at decoupling; 
the amplitude of the polarization is roughly 10% of the pure 
temperature one, as the ordinary Gaussian perturbations. 
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To compute the CMB power spectra induced by a pop- 
ulation of bubbles, we simulate 10° x 10° CMB sky patches 
in which they are placed randomly and fill a fraction x of 
the volume. By taking the angular power spectrum, we ana- 
lyze first the impact of the bubbles only on polarization and 
temperature. We find the signal lying entirely on sub-degree 
angular scales, characterized by a series of narrow peaks for 
polarization and a broader one for temperature, reflecting 
the signals from each bubble singularly. The amplitudes of 
these peaks depends in a rather simple manner by volume 
fraction, dimension and central density contrast of the bub- 
bles, from which we infer an empirical formula giving the 
approximate bubbly CMB power; we also use this formula to 
constrain the bubbly distribution parameters starting from 
existing sub-degree CMB observations. When considering 
bubbles together with Gaussian fluctuations from slow roll 
inflation, as realistical if a bubble nucleating era occurred 
before the end of inflation, we find Gaussian and bubbly 
CMB power spectra approximatively add because of their 
uncorrelation. 

Bubbly density perturbations are remnant of inflation- 
ary tunneling phenomena. Their detection would bring an 
invaluable insight into the hidden sector of high energy 
physics. In this work CMB bubbly signals have been com- 
puted and could be compared with the data from the future 
high resolution observations provided by the MAP and the 
Planck missions in the next decade. 
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Figure 1. Radial behavior of the bubbly temperature perturba- 
tion at different times. The initial condition remains unchanged 
until the horizon reenter occurs. The successive oscillations gen- 
erate CMB perturbation propagating outward with the photon- 
baryon fluid sound velocity. 
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Figure 3. Polarization (up) and temperature (down) anisotropics 
from bubbles with the indicated radii intersecting exactly the last 
scattering surface. The sound waves have been photographed by 
the decoupling photons in both cases. For the polarization, the 
central perturbation spot is absent. Thin long(short) dashed lines 
are anisotropics from bubbles placed at a distance +(— )R with 
respect to the last scattering surface. Thin continue lines in the 
bottom panels show the signals from the two first temperature 
multipolcs 0o and 0i. 



Figure 2. Radial behavior of the bubbly polarization perturba- 
tion at different times. The initial condition remains unchanged 
until the horizon reenter occurs. The successive oscillations gener- 
ate CMB polarization waves propagating outward as in the pure 
temperature case. The substantial difference is the absence of 
central perturbation, because of the spherical symmetry of the 
perturbations. 
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Figure 4. Realistic picture of the bubbly CMB polarization 
anisotropies. The simulated patch is 2° X 2°. The bubble has ra- 
dius R = 20h~ 1 Mpc and central density contrast at decoupling 
8 = 10~ 2 (bottom) and <5 = 5 X 10~ 3 (up). Together with the 
bubbly signal, the portion of sky simulated contains a pure CDM 
Gaussian one. 




Figure 5. Realistic picture of the bubbly CMB temperature 
anisotropies. The simulated patch is 2° X 2°. The bubble has ra- 
dius R, = 20/i -1 Mpc and central density contrast at decoupling 
S = 10~ 2 (bottom) and <5 = 5 X 10~ 3 (up). Together with the 
bubbly signal, the simulated portion of the sky contains a pure 
CDM Gaussian one. Note the presence of the central dark spot, 
absent in the polarization case. 
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Figure 6. The theoretical (dashed line) CMB power spectra for 
the CDM Gaussian perturbations model adopted in this work has 
been simulated and reconstructed (solid line) in a 10° X 10° sky 
patch. Temperature (up) and polarization (down) agree quite well 
for I > 50. 



Figure 7. CMB temperature (up) and polarization (down) power 
spectra for bubbles with R = 20h~ 1 Mpc, S = 10 -3 and filling 
volume fractions between 10% (thin line) and 50% (heavy line). 
Note how the signal is confined to sub-degree angular scales / > 
200. 
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Figure 8. CMB temperature (up) and polarization (down) power 
spectra for Gaussian CDM perturbations and bubbles with R = 
20h~ 1 Mpc, S = 10~ 3 and filling 50% of the volume (heavy line). 
Also plotted are the purely Gaussian perturbations (thin line) 
and the existing experimental data on sub-degree angular scales 
I > 200. 
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Figure 9. CMB temperature (up) and polarization (down) power 
spectra for Gaussian CDM perturbations and bubbles with R = 
20/t- 1 Mpc,<5 = 3 ■ 10~ 3 and filling 10% of the volume (heavy 
line). Also plotted are the purely Gaussian perturbations (thin 
line) and the existing experimental data on sub-degree angular 
scales I > 200. 
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